Thermal correction to resistivity of 2D electron (hole) gas in low-temperature 

measurements at B =0 
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We calculate the zero magnetic field resistivity, taking into 
account the degeneracy of the 2D electron (hole) gas and the 
thermal correction due to the combined Peltier and Seebeck 
effects. The resistivity is found to be universal function of 
temperature, expressed in units of -^(kpl) -1 ■ 

PACS numbers: 73.40. Qv, 71.30+h, 73.20.Fz 



Recently, a great deal of interest has been focussed on 
the anomalous behaviour oLa wde variety of low-density 
2D electron tr U and hole U" u systems, whose resistiv- 
ity unexpectedly decreases as the temperature is lowered, 
exhibiting a behaviour generally associated with metals, 
rather than insulators. In particular, it has been found 
that, below some critical 2D electron density n c s , cooling 
causes an increase in resistivity, whereas at n s > n c s the 
resistivity decreases. Although numerous theories have 
been put forward to account for this effect, the origin 
of this metallic behaviour is still the subject of a heated 
debate. 

In the present paper, we report on a study of low- 
temperature transport in 2D electron gas at zero mag- 
netic field, taking into account both the electron de- 
generacy aHdj-.the Peltier-effect-induced correction to 
resistivity. 13' It is well known that ohmic measure- 
ments are carried out at low current density in order to 
prevent heating. Usually, only the Joule heat is con- 
sidered to be important. In contrast to the Joule heat, 
the Peltier and Thomson effects are linear in current. 
As shown in ||, 0, the Peltier effect influences ohmic 
measurements and results in a correction to a measured 
resistance. When current is flowing, one of the sam- 
ple contacts is heated, and the other cooled, because of 
the Peltier effect. The established temperature gradi- 
ent is proportional to the current. The Thomson heat is 
then proportional to squared current and can therefore 
be neglected. Then, the voltage drop across the circuit 
includes the thermoelectromotive force induced by the 
Peltier effect, which is linear in current. Finally, there 
exists a thermal correction Ap, to the ohmic resistivity, 
p, of the sample. As was demonstrated in 0, for degen- 
erate electrons, Ap/ p m (kT/p) 2 , where p is the Fermi 
energy. Hence, the above correction may be comparable 
with the ohmic resistance of a sample when kT ~ p. 
We further discuss the features of thermal correction 
within 2D electron-density-modulated low-temperature 
ohmic measurements, ur 

Let us consider for clarity a 2DEG sample and dc cur- 
rent flowing in it. The 2DEG structure (MOS, quan- 



tum well, etc.) is arbitrary, electrons are assumed to 
occupy the first quantum-well subband with isotropic 
energy spectrum e(k) — 4^-- Here, m is the electron 
effective mass, p=fik is the electron quasi-momentum, 
and k is the wave vector. The sample is connected (see 
Fig.|f|, insert) by means of two identical leads to the cur- 
rent source. Both contacts are assumed to be ohmic. The 
voltage is measured between the open ends ("c" and "d") 
kept at the temperature of the external thermal reservoir. 
The sample is placed in a sample chamber (not shown) 
with mean temperature To. 

According to our basic assumption, the contacts ( "a" 
and "b") may have different respective temperatures T a 
and Xf,. With the temperature gradient term included, 
the current density j and the energy flux density q are 
given by 



j = cx(E-aVT), 

q= (oT-C/e)j-xVT, 



(1) 



Here, E = VC/e is the electric field, and C, = p — eip 
is the electrochemical potential. Then, a = Ne 2 T/m 
is the conductivity, N is the 2D electron concentration, 
r is the momentum relaxation time, h is the thermal 
conductivity, and a is the 2DEG thermopower. 

It is well known that the Peltier heat is generated by 
current flowing across the interface between two differ- 
ent conductors. At the contact ( "a" in Fig.l], insert ), 
the temperature T a , electrochemical potential £, normal 
components of the current I = jd, and energy flux qd 
are continuous. Here, d is the sample width. Then, there 
exists a difference of thermopowers Aa = a me — a, where 
a me is the thermopower of the metal lead. For Aa > 
and the current direction depicted in Fig.[l], contact "a" 
is heated, and contact "b" is cooled. Thus, the contacts 
are at different temperatures, and T a — Tb = AT > 0. 

In general, one can easily solve Eq.(|jJ), and then find 
AT for an arbitrary circuit cooling. Since the electron- 
phonon coupling is weak below ~1K, the heat conduction 
from 2DEG to mixing chamber could predominately oc- 
cur through the ohmic contacts of the sample and the 
leads connected to them. However, the experimental ob- 
servations in Ref. || demonstrate that the electron gas 
is, in fact, the dominant thermal resistance in this prob- 
lem. Actually, the cooling of 2D electrons with respect to 
bath is provided by thermal conductivity found to follow 
Wiedemann-Franz law. Accordingly, we neglect further 
the contact related cooling of 2D electron gas. Then, for 
actual / — > case we will omit the Joule heating. Let 
us suppose for a moment that the cooling conditions are 
adiabatic, with the 2D electron gas thermally insulated 
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from the environment. These assumptions will be later 
justified for an actual 2DEG system. We emphasize that 
under the above conditions, the sample is not heated. 
Indeed, at small currents, T a sa Tf, fts T . Hence, the 
amount of the Peltier heat, Q a — IAaTo, evolved at con- 
tact "a" and that absorbed at contact "b" are equal. If it 
is recalled that the energy flux is continuous at each con- 
tact, the difference of the contact temperatures is given 
by AT = IAaTolo/xd, where Iq is the sample length. As 
expected, AT is linear in current. 

As shown in Ref. j7j, standard ohmic measurements 
always result in a thermal correction to the resistance 
measured. Using Eq. ([!]), we find for the voltage drop U 
between ends "c" and "d" 



U = RI 



adT, 



(2) 



where R is the total ohmic resistance of the circuit. The 
second term in Eq.(Q) coincides with the expression for 
the conventional thermoelectromotive force, St, under 
zero current conditions. Let us assume that the temper- 
ature gradient is small. In this case <r, a, and k can be 
considered constant. The thermoelectromotive force is 
then given by Et — AaAT. Since AT ~ /, there al- 
ways exists a thermal correction to the ohmic resistance 
AR = e T /I, Finally, the total resistivity of the 2DEG- 
sample is given by LP 



1 



a 
1 



(3) 



where p = \ja is the ohmic resistivity of the sample. 
In Eq.(|J) we take into account that for the actual case 
of metal leads Aa ~ —a. Then, for the low-temperature 
case in question we omit the phonon-related contribution 
to thermal conductivity. Therefore, x = LTa, where 

2.2 

L = is the Lorentz number. It is noteworthy that 
the above result is valid for a 2D hole gas as well. @ 

Using the conventional Gibbs statistics U and the en- 
ergy spectrum specified above, the 2DEG concentration 



N 



yields 



where f2 



exp(^ 



fcT 



(4) 

)) is the thermody- 



-fcT£ln(l 

k 

namic potential of the 2D electron gas, and £ = kT/p is 
dimensionless temperature. Then, Ao ~ ^jj^ is the den- 

oo 

sity of strongly degenerate 2DEG, and Fn(z) = J x n [l + 

o 

exp(x — z)]~ 1 dx is the Fermi integral. In Fig.[l| we plot 
the temperature dependence of the dimensionless con- 
centration n — N/Nq given by Eq.(^). In the classi- 
cal Maxwell-Boltzman limit (£ < 0, [f| < I), the 2D 
electron density is thermally activated, and, therefore, 
n = |£| exp(— 1/ |£|). In the case of strongly degenerate 



electrons (£ <C I), we obtain n = 1 +£exp(— 1/£). Then, 
at elevated temperatures £ >• 1, the dependence of the 
2DEG concentration n = 1/2 + £ In 2 becomes linear in 
temperature. It is noteworthy that the 2D electron con- 
centration A may exceed the zero-temperature value Nq. 

We emphasize that the 2D electron density is a mono- 
tonic function of temperature (see Fig.^J). Therefore, one 
might expect that the ohmic resistivity p(T) ~ 1/A de- 
creases with increasing temperature at constant carrier 
mobility. We now demonstrate that the total resistivity 
specified by Eq.(||) can, nevertheless, increase in a certain 
temperature range owing to the Peltier effect-related cor- 
rection. 

Following the conventional Boltzman equation formal- 
ism, the explicit formulae for the 2DEG thermopower 
(for the 3D case, see Pisarenko, 1940) can be written as 
follows 



2Fi(l/0 1 
Fo(l/0 £ 



(5) 



Here, we assume, for simplicity, that the electron scat- 
tering is characterized by energy-independent momen- 
tum relaxation time. For strongly degenerated 2DEG 
(0 < £ <C 1), we obtain the temperature dependence 
of the thermopower (Figj|,b) as a = -|[7r 2 £/3 - (1 + 
3£) exp(— 1/£)]. At elevated temperatures (£ > 1) the 
thermopower first grows with temperature, and then ap- 
proaches an universal value a s = ~~ 6 ^ n2 ■ I n the classi- 
cal Maxwell-Boltzman limit (£ < 0, [£| < I) the ther- 
mopower is given by the conventional formulae a = 
— -(2 — l/£). Its worth noting that in the two last cases 
the thermopower is of the order of k/e. Accordingly, the 
thermal correction to resistivity Ap = pa 2 / L may be of 
the order of the 2DEG ohmic resistivity. 

In Fig j^, we plot the temperature dependence of the 
2DEG resistivity given by Eq. (^) at different Fermi tem- 
peratures Tf = p/k. At fixed temperature, the resistiv- 
ity decreases with increasing 2DEG degeneracy. Then, 
for a fixed Fermi energy (e.g., Tp — 0.5 in FigJ| ) the T- 
dependence of the resistivity exhibits metallic behaviour 
at T < Tp, and then becomes insulating ( i.e.^ < ) 
at T > Tp. Within the low-temperature metallic region 
£ <C 1, the 2DEG resistivity can be approximated (see 
dashed line in Fig.| ) with p tot = p Q (l + tt 2 £ 2 /3), where 
Po=-^i{kpl)^ 1 is the resistivity at T 
is the Fermi vector, and / 



0, kp = yJ2mp/h 
hkpT/m is the mean free 
path. Then, for the high-temperature (£ > 1) insulating 

-Po |Tirfw72 > shown 



region we obtain the asymptote p*" 



m Fig.|bydotteoUi ne. This result is confirmed by recent 
experiments 113" G-l shown that for temperatures well be- 
low the Fermi temperature the metallic region data obey 
a scaling law where the disordered parameter jkpl appears 
explicitly These experimental observations O therefore 
rule out interactions, the shape of the potential well spin- 
orbit effects as possible origins of the metallic behaviour 
mechanism. Then, according to Ref. |0| the "metallic" 
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state can not be associated with e-e induced quantum in- 
terference effects. We argue that the semi-classical mech- 
anism discussed above may be responsible for observed 
T-behaviour of 2D resistivity. 

Let us analyze in more detail the cooling conditions, 
which are known to influence the thermal correction to 
resistivity. O' It will be recalled that in the case of 
adiabatic cooling the electron temperature differs from 
the bath temperature T . We now consider the op- 
posite situation of electron cooling due to, for exam- 
ple, finite strength of electron-phonon coupling. Let the 
phonons are thermalized with respect to bath tempera- 
ture T . Following Ref. g5|, below -0.6K in Si MOS- 
FET's the electron-to-phonon thermal exchange is given 
by a(T 3 - T 3 ), where a =2.2 x I(T 8 W/K 3 cm 2 . When 
T — Tq <C To, the thermal correction to the resistiv- 
ity Ap is suppressed by the factor (3 = tal ^ h x , where 

A = loToy^3a/A>c is a dimensionless parameter. Actually, 
A is the ratio of outgoing and internal heat fluxes asso- 
ciated with phonon related thermal leakage and electron 
heat diffusion, respectively. 

When A « 1, the local cooling due to phonons can 
be neglected and the adiabatic approach is well justi- 
fied. In the opposite case of intensive cooling (A ^> f), 
the difference of the contact temperatures AT becomes 
smaller, and, therefore, the thermal correction to re- 
sistivity Ap vanishes. For Iq = Imm, To = 50mK, 
a = 2e 2 /h = 8 x 10 _5 Ohm _1 ( typical critical region 
conductance) we obtain A = 1.2, and, therefore, (3 = 0.8. 
It worthwhile to notice that Peltier correction to resistiv- 
ity becomes greater at ultra-low temperatures for short 
samples, since A ~ ^oT). In real experiments both the 
electron-phonon coupling and the sample-to-bath ther- 
mal exchange may be important. 

We emphasize that both dc and ac ohmic measure- 
ments lead to a thermal correction. The correction is, 
however, strongly damped at high frequencies because of 
the thermal inertial effects. As demonstrated in Ref. ||, 
the above quasi-static approach is valid below some crit- 
ical frequency f cr — xAo- For example, for degenerate 
electrons the thermal diffusion coefficient \ is of the order 

of the diffusion coefficient D = ^ {^j^j ■ Assuming 

a = e 2 /h, l = 1mm, for GaAs-based structure we ob- 
tain x ~ h/m, hence f cr — 1. 5kHz. We suggest that the 
spectral dependence of the 2D resistivity can be used to 
estimate the thermal correction. 

In conclusion, low-temperature ohmic measurements 
of a 2D electron (hole) gas at B = involve a thermal 
correction caused by the Peltier effect. The magnitude 
of thermal correction depends on the 2DEG degeneracy 
and actual cooling conditions. The resistivity of 2DEG 
with thermal correction included is found to be universal 
function of temperature, expressed in units h/e 2 (kpl)~' 1 . 
This universal behaviour correlates with that found in 
experiments. 
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FIG. 1. Themperature dependence of 2DEG concentration 
(a) and thermopower (b) given by Eq.(^) and Eq.(^) respec- 
tively. Asymptotes: |£| -C 1 - dotted line, £ S> 1 - dashed line. 
Inserts: the experimental setup (top); position of the Fermi 
level with respect to the bottom of quantum-well subband 
(bottom). 



FIG. 2. T-dependence of the 2DEG resistivity, given by 
Eq.(|H§ for T F = 2; 1.75; 1.5; 1.25; 1; 0.75; 0.5; 
0.3; 0.25; 0.2; 0.15; 0.1; 0.05; 0.01K. Asymptotes: £ < 1 - 
dashed line, £ > 1 - dotted line for Tf = 0.5K. Insert: density 
dependence of the 2DEG resistivity within the T =0.5-0.9K 
range. 
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